Multistability, noise and attractor-hopping: 
The crucial role of chaotic saddles 
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We investigate the hopping dynamics between different attractors in a multistable system under 
the influence of noise. Using symbolic dynamics we flnd a sudden increase of dynamical entropies, 
when a system parameter is varied. This effect is explained by a novel bifurcation involving two 
chaotic saddles. We also demonstrate that the transient lifetimes on the saddle obey a scaling law 
in analogy to crisis. 
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Systems with a large number of coexisting stable states 
have been the subject of increasing interest recently. This 
multistable behavior occurs in many different fields like 
optics chemistry , neuroscience 0] , semiconductor 
physics^ , plasma physics [1| and coupled oscillators (Qj . 
Moreover, addition of noise to multistable systems has 
led to several interesting phenomena, like noise-induced 
preference of attractors , directed diffusion , noise- 
enhanced multistability and chaotic itinerancy 
The latter effect, which consists of hopping between dif- 
ferent attractors caused by the noise, has also been ob- 
served experimentally in an optical system [l^ . 
In this Letter we focus on the dynamics of the attractor- 
hopping process by taking a simple model as a paradigm. 
Due to the noise the attractors become 'metastable' and 
the trajectory starts hopping between the different at- 
tractors. In the case of fractal basin boundaries, this 
hopping process consists of two steps. In the first one 
the trajectory leaves the open neighborhood about the 
attractor according to Arrhenius law (13i |. in the second 
one, the trajectory bounces around on the chaotic saddle 
before cascading again into the open neighborhood of an 
attractor. The general mechanism for the first step of 
the hopping process is the same for bi- and multistable 
systems. In both cases we find Arrhenius law; quanita- 
tive differences are due to different relative sizes of basins 
of attraction However, for the attractor-hopping dy- 
namics, our study points out a major difference between 
bistable and multistable systems: While in bistable sys- 
tems the structure of the saddle separating the two at- 
tractors (saddle point or chaotic saddle) does not play 
a role for the hopping characteristics, it is essential for 
the attractor-hopping in a system possessing a multitude 
of attractors. The structure of the chaotic saddles influ- 
ences the nature of the hopping process, in particular it 
determines which transitions between attractors are pos- 
sible. Bifurcations in the chaotic saddles lead to changes 
in the 'accessibility' of attractors in the hopping process. 
Employing symbolic dynamics, we assign one symbol to 
each attractor, thus transforming our hopping time series 
onto a symbolic string. As a next step we compute both, 



the Shannon and the topological entropy, which can be 
regarded as measures of complexity for the hopping dy- 
namics. As a system parameter is varied the complexity 
of the hopping process changes beyond a certain thresh- 
old value. We show, that this change can be related to a 
novel bifurcation, namely a merging of two chaotic sad- 
dles accompanied by the emergence of additional points 
filling the gap between the formely separated saddles. 
This bifurcation is mediated by a snapback repellor |3| i 
whose eigenvalues determine the scaling law for the tran- 
sient lifetimes on one chaotic saddle close to the bifurca- 
tion point. 

Our basic prototype model of multistablity, which cap- 
tures the main features of highly multistable systems, 
namely more than two final states (attractors) and a 
complexly interwoven fractal basin boundary separating 
these states, is given by the 10-fold iterate of two coupled 
logistic maps: 
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We fix 7, the coupling strength, at 7 = 0.29, and vary a, 
the nonlinearity, in the range a G [0.72,0.755]. For this 
parameter set, there exists a stable period 10-orbit for 
the two coupled logistic maps themselves. Since we con- 
sider the 10-fold iterate, our map exhibits 10 coexisting 
fixed point attractors, where each 5 lie above and below 
the symmetry axis x — y, respectively. Using only ini- 
tial conditions below the symmetry axis, we can restrict 
our study to a system possessing 5 coexisting fixed point 
attractors, 5 repellors with two unstable directions and 
10 saddle points. Each of the 5 attractors is surrounded 
by an open neighborhood of different size, in which all 
initial conditions converge to the corresponding attrac- 
tor inside this neighborhood. The basins of attraction 
as a whole have a complex fractal structure due to a 
homoclinic bifurcation which occurs already for smaller 
values of the nonlinearity a. As a next step we apply 
noise to the system, using Gaussian, white noise with 
standarddeviation a added to the x and y component of 
the 10-fold iterate of Q). Small noise causes the system 
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now to alternate between the different states, where long 
periods close to a fixed point, comparable to laminar mo- 
tion, are interrupted by short, sudden bursts, which are 
reminiscent to intermittent behavior. During the bursts 
the motion takes place on the chaotic saddles separat- 
ing the attractors, until the trajectory is again being in- 
jected in the neighborhood of one of the 5 fixed points. 
This complex dynamics is depicted in Fig. ^ According 
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FIG. 1: Noisy time series of equation Q for a = 0.73,7 — 
0.29 and a — 0.012. There are clearly 5 distinct almost pe- 
riodic states, marked with roman numbers on the right hand 
side. These states are interrupted by bursts, where the motion 
takes place on a chaotic saddle. 



to our aim we are now interested in the characteristic 
properties of the hopping dynamics. Thus we address 
the question, how the complexity of the hopping dynam- 
ics changes with the nonlinearity parameter a, namely 
whether every fixed point has a positive transition prob- 
ability to every other fixed point, via a transient on a 
chaotic saddle. For this purpose, we employ the concept 
of symbolic dynamics. We use an encoding scheme in 
which we assign one symbol to each attractor Ne- 
glecting the number of iterations the trajectory spends 
close to a fixed point, a symbol is given for every fixed 
point and only after a jump out of an attractor a new 
symbol is bestowed, according to the attractor where the 
trajectory lands at. Using this scheme, we focus only on 
the structural properties of the jumps, not taking fully 
into account the complete temporal evolution in each it- 
eration. As a quantitative measure of the complexity of 
the symbol string we use the Shannon entropy, in analogy 
to the Kolmogorov-Sinai entropy |1(tI | , given by 
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of occurrence and ff„ the block entropy of block length n. 
Numerically, the quantity lim„^oo(-ffri-i-i — ^^n) converges 
already for n = 1, indicating that the hopping dynamics 
is a Markov process of first order. This corresponds to the 
intuitive expectation, that the dynamics possesses only 
memory of the last state it was dwelling on. Higher cor- 
relations are suppressed by the long laminar-like motion. 
Consequently a transition matrix between the different 
fixed points can be constructed, where the entries are 
the probabilities for a transition between two states. 
Another interesting quantity to investigate is the topo- 
logical entropy. One way of computing it is to built the 
Stefan transition matrix [l^ . where the entries of the 
5x5 matrix are 1 and 0, depending on whether a tran- 
sition took place or not, respectively. The logarithm of 
the largest eigenvalue of this matrix yields the topolog- 
ical entropy. In our numerical implementation, we use 
a certain, very small cutoff limit of 0.001 for the tran- 
sition probabilities, below which we regard the value as 
zero. Thus we neglect transition, having extremely small 
probabilities compared to the other transitions. How- 
ever, the exact value of this cutoff does not change the 
results significantly, as long as it is small enough. The 
entropies shown in Fig. El remain unchanged for cutoff 
limits between 10"'^ to 10~^. It is important to note 
that, although the escape times out of the stable states 
depend on the noise level, the transition probabilities re- 
main constant for a wide range of noise values. Thus the 
considered entropies are not affected by the noise ampli- 
tude for a rather large interval of noise levels. In Fig.Elwe 
present the evolution of these two quantities as the non- 
linearity a changes. The curve for the topological entropy 
appears to be a monotonic curve with a devil's staircase 
like behavior Every time it increases, at least one 




where S = SiS2---Sn denotes a finite symbol sequence con- 
sisting of n elements Si = 1,2,..., 5, p{S) its probability 
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FIG. 2: The Shannon entropy (o) and the topological entropy 
(□) of the symbol sequence generated from equation Q under 
variation of the nonlinearity parameter a. Both quantities are 
normalized to 1. 



new transition between hitherto unconnected fixed points 
is created. Below the homoclinic bifurcation a < 0.72, 
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where the basin boundaries are smooth, only one tran- 
sition per fixed point is possible, as the fixed points are 
located on a closed curve made up by the unstable man- 
ifolds of the saddlepoints. These unstable manifolds im- 
pose a direction. Consequently, the topological entropy is 
zero. Above a value of a « 0.75, all transitions can occur 
and the maximum of the topological entropy is reached. 
The Shannon entropy is always less than the topological 
entropy. Besides the overall increase of both entropies, 
the plateau between 0.74 < a < 0.744 suggests, that 
there is a change in the dynamical behavior beyond it. 
This change is due to a change in the topological struc- 
ture of the system and relates to a new bifurcation as we 
show next. 

We compute the chaotic saddle, responsible for the 
chaotic dynamics, by the PIM-triple (proper interior 
maximum) method [l9j. The result is depicted in Fig. 
13 and 21 where the chaotic saddles are shown for two 
different values of a together with the 4-fold preimages 
of one repellor R marked with a □. Slightly above 
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FIG. 3: Two chaotic saddle rings in different gray scale to- 
gether with the 4-fold preimages of the repellor K marked 
with a □ for a = 0.740. 



the homoclinic bifurcation (a = 0.725) the saddles are 
very thin, consisting of two separate rings. At a value of 
OL ~ 0.755, there exists only one large piece of the nonat- 
tracting chaotic set (not shown). This results from the 
saddle merging, which took place at a w 0.7430 (Fig.^. 
The two bands of chaotic saddles merge and the repellor, 
formerly located ouside of the saddles, is now embedded 
in the large saddle, including some of its preimages. The 
merging is accompanied by the successive filling of the 
formerly empty gap between the two saddles. A similar 
effect has been observed in p3| where it is called 'spilling', 
in chaotic scattering [ij , and in the gap filling crisis [22| ■ 
In contrast to the latter this happens smoothly in the 
saddle merging. 

However, the actual bifurcation takes place at a slightly 
lower value of a w 0.7428 (not shown). At this param- 
eter value some of the preimages of the marked repellor 
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FIG. 4: Chaotic saddle together with the 4-fold preimages 
(-I-) of the repellor R marked with a □ for a = 0.7430. 



touch its unstable manifold, which constitutes the border 
of the chaotic saddle. Through this mechanism, a snap- 
back repellor has developed. A fixed point p is called a 
snapback repellor, if i) all eigenvalues of p have absolute 
values larger than 1, ii) there exists a, q G Wl^^{p), the 
unstable manifold of p, such that the M-fold iterate of 
the map {q) =_p for some positive integer M and iii) 
det DJ^^^{q) 7^ 0[lj|. While conditions i) and iii) are true 
for all a considered here, it is condition ii) that becomes 
fulfilled during the bifurcation, as the unstable manifold 
of the repellor touches its preimages. This bifurcation 
is responsible for the emergence of the additional points 
in the gaps, (see Fig. At a shghtly higher value of 
a fa 0.7430 the preimages touch not only the unstable 
manifold, constituting the border of the chaotic saddle, 
but the saddle itself. As the chaotic saddle is an invari- 
ant set, the repellor becomes thus embedded in the saddle 
and the merging of the two rings has taken place. This 
merging has an important consequence for the dynam- 
ics of the attractor-hopping. Not only jumps between 
neighbors, but arbitrary jumps become, in principle, pos- 
sible. Above the critical value, the connected pieces of 
the chaotic saddle enable the trajectory to jump between 
the formerly separated parts, thereby causing a sharp in- 
crease in the number of the possible transitions, as can 
be seen from the topological entropy above a — 0.7430 
in Fig. El 

Since the bifurcation described above can be considered 
as a kind of crisis, it is interesting to see, whether the 
scaling laws for the characteristic transient lifetimes |2^ 
apply in this case as well. We compute the average es- 
cape times out of one (still separated) chaotic saddle ring 
slightly below the bifurcation. In Fig.jSlthe data are plot- 
ted, resulting in a scaling relation: 

(r)^ (a -«,)-% (3) 

where etc = 0.7428. The best fit to the data yields an 
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exponent e = 1.86. That value can be reproduced by 

e = l/2 + (ln(/3i)/ln(/32)) (4) 
e = 1/2 + (ln(2.45)/ln(1.9)) = 1.85, (5) 

where Pi and (32 are the expanding eigenvalues of the 
mediating repellor. This coincidence has been con- 
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FIG. 5; Scaling of the transient lifetime to stay in one ring of 
the chaotic saddle with the nonlinearity a and Qc = 0.7428. 
The slope of the log-log plot is e = 1.86. 



firmed also for other values of the coupling strength 
7 = 0.285, 0.287, 0.292. It is worth mentioning that in the 
original derivation of this formula, which applies to the 
heteroclinic tangency crisis. Pi and P2 are the expanding 
and contracting eigenvalues of the mediating saddle point 
[i^ l. However, in our case both. Pi and P2, are expanding 
eigenvalues of the mediating repellor. We conjecture this 
relation to hold for the general case of a crisis caused by 
a snapback repellor. The reason, that a law for the het- 
eroclinic tangency is supposed to hold in a homoclinic 
case of a snapback repellor is, that the mechanism for 
approaching the fixed point is different here. The fixed 
point is not reached by successive iterations, but sud- 
denly with a jump of noninfinitesimal size from a finite 
distance. 

In conclusion, we have studied a simple paradigmatic 
map possessing a multitude of coexisting stable states un- 
der the influence of small noise. We have shown that the 
nature of the attractor-hopping process depends, in con- 
trast to the bistable case, crucially on the structure of the 
chaotic saddles separating the attractors. Bifurcations of 
the saddles lead to a change in the hopping dynamics 
which is manifested in a sudden increase of dynamical 
entropies. In particular we described a novel bifurcation, 
a merging of two chaotic saddles involving a snapback 
repellor, which results in a change of the 'accessibility' 
of the attractors. Finally, a conjecture for the scaling 
law of the transient lifetimes on the chaotic saddle has 
been provided. The hopping dynamics, the bifurcation 
and the scaling law should be observable in experimental 
realizations of such systems. We expect this bifurcation 



also to occur in systems with three or more dimensions, 
where the mediating saddle has at least two unstable di- 
rections. 
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and the DFG for financial support. 
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